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Abstract 

We have magnetically imaged interlayer Josephson vortices emerging from 
an ac face of single crystals of the single layer cuprate high-T c superconduc- 
tor (Hg,Cu)Ba2Cu04 +< 5. These images provide a direct measurement of the 
c-axis penetration depth, A c ~ 10 /im. This length is a factor of 10 longer 
than predicted by the interlayer tunneling model for the mechanism of super- 
conductivity in layered compounds, indicating that the condensation energy 
available through this mechanism is 100 times smaller than is required for 
sup er conductivity. 
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In the interlayer tunneling (ILT) model for superconductivity in layered superconductors 
such as the cuprates []]]-§]], transport of carriers between the planes is incoherent in the 
normal state, but coherent interlayer transport is allowed for Cooper pairs. The coherent 
pair tunneling lowers the c-axis kinetic energy, supplying the superconducting condensation 
energy, E c . Anderson ||||, Leggett H0], and Chakravarty [[J have each argued that the 
comparison of the experimentally measured c-axis penetration depth, A c , and the value 
determined within the ILT model from the condensation energy, Xilt, is an important test 
of the ILT mechanism. In all presently published versions of the theory, A c ~ Xilt- 

The calculation of Xilt and its comparison with experiment are most straightforward in 
materials with a single copper oxide layer per unit cell, such as La 2 _ x Sr x Cu04 +( 5 (Xilt — 3/im 
§), HgBa 2 Cu0 4+5 (Xilt ^ l^m §), and Tl 2 Ba 2 Cu0 6+5 (Xilt ^ l/^m §). Different 
experiments on these materials have not provided a single answer to the question 0,0. In 
La 2 _ 3; Sr 3 ;Cu04 +( 5 (La-214), measurements of the Josephson plasma frequency lo p = cA c 7 1 e -1 / 2 
(where c is the speed of light and e is the dielectric constant of the interlayer medium) 
PJTT|] are in good agreement with the predictions of the ILT model ||. In Tl 2 Ba 2 Cu0 6+< 5 
(Tl-2201), Moler et al. observed interlayer Josephson vortices with A c = 17 — 21 microns 
[TIJ. Measurements of the plasma resonance u p and dielectric constant e by Tsvetkov et al. 



indicate A c = 17 microns In HgBa 2 Cu04 + 5 (Hg-1201), Panagopoulos et al. obtained 

a value of A C (T = 0) = 1.36±0.16/zm from magnetic susceptibility data on oriented powders 
13| . As both Anderson and Leggett have pointed out, a confirmation of the c-axis 



penetration depth in Hg-1201 seems essential to resolving this issue. 

In this Letter we directly measure the c-axis penetration depth in Hg-1201 by magnet- 
ically imaging interlayer Josephson vortices emerging from the ac face of single crystals. 
These measurements give A c ~ 10/im, much longer than the value A c ~ lfim reported by 
Panagopoulos et al. fD3f . Our value for A c , when combined with previous results from Tl- 
2201, make it appear unlikely that presently published versions of the ILT model are the 
correct description for superconductivity in the cuprate high-T c superconductors. 

The Hg-1201 crystal growth has been described previously [14]]. Our method consists 



in preparing a Ba/Cu/O precursor in flowing oxygen, and then mixing it with HgO to 
make a 0.8:2:1.2 ratio of Hg:Ba:Cu. An alumina crucible containing this powder is then 
sealed in a silica tube. After 48 hours of thermal treatment, black, platelet-like crystals 
are extracted. The largest crystals (typical dimensions ~lx 1x0. 08mm) are selected for 
transport and SQUID imaging measurements. EDX and electron diffraction measurements 
lead to the formula Hgo.sCuo.2 for the mixed mercury layer. The copper position is displaced 
with respect to the Hg site in the mixed layer as shown from structural refinements. No 
intergrowth nor extended defects have been shown by high resolution microscopy [|lf|]. The 
Hg-1201 single crystals were mounted in epoxy so that the ac face was aligned vertically, 
and then the sample and epoxy were polished to make a flat, smooth surface to scan. 

The magnetic imaging measurements were made with a scanning SQUID microscope 
|T6| , in which a sample is scanned relative to a superconducting pickup loop oriented nearly 
parallel to the sample surface. The data are represented as the magnetic flux $ s , the integral 
of the normal component of the magnetic field over the pickup loop, vs the position of the 
pickup loop in the x, y plane. The pickup loop is fabricated with well-shielded leads to 
an integrated Nb-Al 2 3 -Nb SQUID. For these measurements, we used a SQUID with our 
smallest pickup loop, an octagon L = 4/xm in diameter with 0.8/im linewidth. The silicon 
substrate carrying the SQUID was polished to a sharp tip a few microns from the edge of the 
pickup loop, the substrate was oriented an angle 9 ~ 20° from parallel to the sample, and 
scanned with the tip in direct contact with the sample. Fits to data on Abrikosov vortices 
indicate that the pickup loop is about ZQ=lixm. above the surface of the sample under the 
experimental conditions used here. Our images were made with both sample and SQUID 
immersed in liquid helium at 4.2K. 

Figure 2 shows a SSM image of a 512/imx256/im area of the crystal and surrounding 
epoxy. The outer edges of the crystal face are indicated by dashed lines. The sample was 
cooled and imaged in a field of about lOmG. Visible in this image are about 50 interlayer 
Josephson vortices. These vortices are remarkably uniform in shape. Some are sufficiently 
isolated from their neighbors that overlapping fields do not interfere with detailed modelling. 



Although one must consider that the vortices may be pinned in areas with unusually weakly 
coupled planes, the uniformity of the vortex shapes over a large area of the crystal face makes 
us believe that we are measuring an intrinsic penetration depth. This belief is supported by 
the good agreement between bulk plasma resonance measurements of uij and local interlayer 
Josephson vortex imaging measurements of A c in Tl-2201 W%. Three vortices chosen for 



further analysis are shown in figure 2. Cross-sections through the image data parallel to the 
layers are displayed in Figure 3. 

The decay of the observed vortex magnetic flux perpendicular to the layers (along the 
c-axis) is determined by the size of the pickup loop. In contrast, the extent of the images of 
the vortices along the layers provides a direct measure of A c . 

For quantitative modelling of the shape of the vortex we use the results of Clem and 



Coffey [17]. Neglecting the influence of the surface on the fields, the z-component of the 



magnetic field of an interlayer vortex is given by |17| : 



B z {x, y,z — 0) — —^—K (R), (1) 

where \ a b is the in-plane penetration depth, K is a modified Bessel function of the second 
kind of order 0, $o — hc/2e is the superconducting flux quantum, h is Planck's constant, 
e is the charge on the electron, R = ((s/2\ a b) 2 + (x/\ a b) 2 + (y/^c) 2 ) 1 ^ 2 , s is the interplane 
spacing, x is the distance perpendicular to the planes, and y is the distance parallel to the 
planes. Since for our experiments s << \ a b << (L,Zq) « A c , we neglect both s and \ a b- 



With this assumption, the fields from Eq. 1 are propagated to a height z = z [[Tj|] and then 
summed over the geometry of the pickup loop. This model has two free parameters: A c , 
which determines the length of the vortex, and Zo, which determines the magnetic amplitude 
of the vortex image. Fits to the three cross- sect ions (Figure 3) yield consistent values for 
the interlayer penetration depth A c = 10 ± \\xm. 

This measured value is about ten times longer than the theoretical value, Xilt = 1 ± 
0.5/zm [^]. Since the ILT condensation energy is proportional to l/\ 2 , the ILT supplied 
condensation energy is about 100 times smaller than the estimated actual energy in Hg- 



1201 §. 

A more conventional estimate of A c comes from the Lawrence-Doniach model [|19|]. For 
diffusive pair transfer (parallel momentum not conserved) between superconducting layers 
p~9^pl| , the Josephson current between two identical superconducting sheets at T=0 is given 
by 12J 

where A(0) is the zero temperature energy gap and i? Cjn is the normal state c-axis inter- 
plane resistance. The interlayer penetration depth is given by Aj_ = (c^q/S^sJq) 1 ^ 2 [|T7 



The application of this model requires at least two assumptions: that R n is temperature 
independent below T c , or at least that the temperature dependence can be understood well 
enough for extrapolation ||, and that the gap is s-wave p2[-f2"4|j . d-wave superconductivity 
would tend to increase A c from this value: in a purely tetragonal o?-wave superconductor, 
with purely diffusive pair transfer, the coupling would be reduced to zero [p2}-|2~4]|. With 
these considerations in mind, one would not expect the Ambegaoker-Baratoff model to ap- 
ply quantitatively, but it is nevertheless natural to look for a correlation between A c and 

Rc,n- 

Figure 4a shows measurements of the c-axis resistivity of two of our Hg-1201 crystals as 
a function of temperature. These measurements were made by evaporating four Ag stripes, 



two on each ab face of the crystal, as diagrammed in the inset of Fig 4a | 25fl . The contacts 
as prepared have high resistance, but after an annealing step at 400°C for 10 minutes the 
contact resistances drop to 1 to 2 Ohms. We have carefully checked that such annealing 
does not alter the superconducting transition by measuring the susceptibility of control 
crystals before and after such an annealing step. Because of the relatively high anisotropy 
of Hg-1201 and large thickness of our crystals (~ 80fim), these crystals are electrically 
thick |[26||. Therefore difficulties in separating the in- and out-of-plane conductivities such 



as reported by Hussey et al. [f27| for YBa 2 Cu 4 08 should not occur. The midpoint T C =94K 
of the p c {T) curves is consistent with the onset T C =95-97K observed magnetically on such 



crystals and confirms that the doping is close to optimal. Taking p cn =cr c ^ = 0.5 Ohm- 
cm, and a BCS gap value 2A=3.54k B T c = 28 meV, leads to a predicted penetration depth 
A c = (hc 2 /ATT 2 Aa c ^ 2 = 8/xm. 



Basov et al. [£8] have previously noted a correlation between the values of the penetration 
depth A c and the far-infrared c-axis conductivity, a c ,FiR- In Figure 4b add our results for 
Hg-1201 and Tl-2201 flTOl , replacing g c ,fir with a CjU just above T c . The dashed line is Eq. 
(2), assuming 2 A = 28meV. This simple model, although it assumes the same A for all 
cuprates, and has the shortcomings outlined above, agrees with all of the measurements to 
within a factor of three. 

It is not clear why our results for A c are nearly a factor of 7 larger than those of 
Panagopoulos et al. There are possible systematic errors in both measurements. In the 
analysis of the imaging experiments, the spreading of the vortex near the surface is neglected. 
This neglect is justified within our experimental resolution on the basis of a consideration of 
the free energy associated with vortex spreading in highly anisotropic superconductors []29 



as well as by the quantitative experimental agreement between A c and lv p in Tl-2201 
For the powder magnetization results, the analysis depends on a very accurate characteriza- 
tion of the distribution of particle sizes; it is possible that this analysis could be skewed by 
the presence of a slightly higher number of large particles than used in the modelling |30| . 



A more likely explanation is that the different samples do in fact have different A c values. 
First, the doping in the two samples is not identical. It has been shown that the transport 
anisotropy of several materials can depend strongly on doping |31| , but it has also been 



shown that the anisotropy in Hg-1201 is nearly constant near optimal doping [j32| . Another 



possibility is the influence of the excess copper on the mercury layer. Neutron studies on 
ceramics f33| , |34]j conclude they also have mixed Hg/Cu layers. It is not clear whether the 



ceramics of Panagopoulos et al. have similar substitutions on the mercury layer. 

If the Hg-1201 crystals used in this study and the Tl-2201 crystals used in the previous 
study are true single-layer materials, then the long c-axis penetration depths A c ~ 10 microns 
for Hg-1201 and A c « 20 microns for Tl-2201 pose a serious challenge to the ILT model. 



Although one can hypothesize structural explanations, perhaps based on the excess copper, 
to make these results consistent with the ILT model, such explanations must be found for 
both (Hg,Cu)Ba2Cu04 + 5 and Tl 2 Ba2Cu06+<5. 

These two results fall into a previously unfilled range of penetration depths in the Basov 
correlation. This correlation, A c (he 2 / An 2 &.a c , n ) l l 2 , works surprisingly well given the 
simplicity of the theoretical model and considering that details of the gap anisotropy, the 
band structure, the scattering, and the temperature dependence of the normal-state c-axis 
resistivity are all unknown in these particular materials. 

We would like to acknowledge M.B. Ketchen for the design, and M. Bhushan for the 
fabrication, of the SQUIDs used in our microscope. We would also like to acknowledge 
useful conversations with PW. Anderson, S. Chakravarty, A.J. Leggett, C. Panagopoulos, 
and D. van der Marel. 



7 



REFERENCES 

[1] J. Wheatley, T. Hsu, P.W. Anderson, Nature 333, 121(1968). 

[2] P.W. Anderson, Physica C 185, 11 (1991); Phys. Rev. Lett. 67, 660 (1991); Science 
256, 1526 (1992). 

[3] S. Chakravarty, A. Sudb0, P.W. Anderson, S. Strong, Science 261, 331 (1993). 
[4] P.W. Anderson, ibid. 268, 1154 (1995). 
[5] P.W. Anderson Science 279 1196(1998). 
[6] A.J. Leggett Science 274, 587 (1996). 
[7] A.J. Leggett, Science 279, 1157 (1998). 



S. Chakravarty, cond-mat 980 102$ , Eur. Phys. J. B, in press. There is a factor of 2 



difference between the estimates of Chakravarty and Anderson for Xht- 
[9] S. Uchida, K. Tamasaku, and S. Tajima, Phys. Rev. B 53, 14558(1996). 
[10] K. A. Moler, et al, ibid 279, 1193 (1998). 
[11] J. Schutzmann et al., Phys. Rev. B 55, 11118(1997). 
[12] A. Tsvetkov et al. to be published. 

[13] C. Panagopoulos et al, Phys. Rev. Lett. 79, 2320 (1997). 
[14] D. Pelloquin, et al., Physica C, 205 (1997). 
[15] D. Pelloquin, M. Hervieu, and B. Raveau, to be published. 
[16] J.R. Kirtley et al. Appl. Phys. Lett. 66, 1138 (1995). 

[17] John R. Clem and Mark W. Coffey, Phys. Rev. B42, 6209 (1990); John R. Clem, Mark 
W. Coffey, and Zhidong Haro, ibid 44, 2732 (1991); A. Gurevich, M. Benkraouda, and 
J.R. Clem, ibid 54, 13196 (1996). 

8 



[18] J.R. Kirtley et al. Phys. Rev. Lett. 76, 1336 (1996). 

[19] W.E. Lawrence and S. Doniach, in Proceedings of the 12th International Conference on 
Low Temperature Physics edited by E. Kanda (Academic Press, Kyoto, 1971), p. 361. 

[20] L.N. Bulaevski, Sov. Phys. JETP 37, 1133 (1973). 

[21] V. Ambegaokar and A. Baratoff, Phys. Rev. Lett. 10, 486 (1963); erratum 11, 104 
(1963). 

[22] M.J. Graf, D. Rainer, and J.A. Sauls Phys. Rev. B 47, 12089 (1993). 
[23] A.G. Rojo and K. Levin Phys. Rev. B 48, 16861 (1993). 

[24] RJ. Hirschfeld, S.M. Quinlan, and D.J. Scalapino, Phys. Rev. B 55, 12742(1997). 

[25] V. Hardy et al. Phys. Rev. B 56, 130 (1997). 

[26] H.C. Montgomery J. Appl. Phys. 42, 2971 (1970). 

[27] N.E. Hussey, et al, Phys. Rev. B 56 R11423 (1997). 

[28] D.N. Basov, et al, Phys. Rev. B 50, 3511 (1994). 

[29] D. Huse, private communication. 

[30] C. Panagopoulos, private communcation. 

[31] for example, Y. Nakamura and S. Uchida, Phys. Rev. B 47, 8369 (1993); Y. Kotaka et 
al, Physica C, 235-240, 1529 (1994); R. Yoshizaki et al, J. Low Temp. Phys. 105, 927 
(1996). 

[32] J. Hofer, to appear in J. Phys. C. 

[33] J.L. Wagner et al, Physica C210, 447 (1993). 

[34] A. Asab, et al, Physica C 255, 180 (1995). 



9 



FIGURES 

FIG. 1. Scanning SQUID microscope image of a 256x512 micron area of the edge of a single 
crystal of Hg-1201, cooled in a field of about lOmG, and imaged at 4.2K. The false color lookup 
table corresponds to a total variation in flux through the SQUID pickup loop of O.SS^o- The 
dashed lines indicate the top and bottom ab faces of the crystal. The box indicates the area that 
is expanded in Figure 2. 

FIG. 2. Expanded view of a 54x100 micron area of the image of Figure 1. The dashed lines 
indicate the paths of cross-sections through the data parallel to the planes displayed in Figure 3. 

FIG. 3. The symbols are cross-sectional data through the 3 representative vortices indicated 
in Figure 2, offset vertically for clarity. The solid lines are fits to the data, as described in the 
text. The fit penetration depths A c are as labelled in the figure, and the effective heights are zq = 
0.6±0.1^ m(A), 0.8±0.1/x m(B), and 0.8±0.2^m(C). The error bars are assigned using a doubling 
of the best-fit x 2 as a criterion. 

FIG. 4. (a) c-axis resistivity vs temperature for two of our crystals. The contact geometry used 
is diagrammed in an insert, (b) Correlation plot between c-axis conductivity and c-axis penetration 
depth, following Basov et al. (Ref. 29). We have included our present data for Hg-1201 and previous 
data on Tl-1201. The dashed line is the prediction for diffuse tunneling between superconducting 
layers, assuming a gap value 2A=28meV. In this figure 124=YBa2Cu40s, 2212=Bi2Sr2CaCu20s, 
p S 2=ps3=Pb2Sr 2 RCu30 8 , 214=Lai. 8 4Sro.i6Cu04,123=YBa 2 Cu 3 O a .. 
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